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Abstract 

We calculate the hydrodynamic time scales for a spherical ultra-relativistic 
shell that is decelerated by the ISM and discuss the possible relations 
between these time scales and the observed temporal structure in 7-ray 
bursts. We suggest that the bursts' duration is related to the deceleration 
time, the variability is related to the ISM inhomogeneities and precursors 
are related to internal shocks within the shell. Good agreement can be 
achieved for these quantities with reasonable, not fined tuned, astrophysi- 
cal parameters. The difference between Newtonian and relativistic reverse 
shocks may lead to the observed bimodal distribution of bursts' durations. 

Subject heading: gamma-rays: bursts- hydrodynamics-relativity 
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1. Introduction 

Gamma-ray bursts (GRBs) are most likely generated during deceleration of ultra- 
relativistic particles. A cosmological compact source that emits the energy re- 
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quired for a GRB cannot generate the observed non thermal burst. Instead it 
will create an opaque fireball (Goodman, 1986; Paczyhski, 1986; Piran, 1994). 
If even a small amount of baryonic matter is present then the ultimate result of 
this fireball will be a shell of ultra-relativistic particles (Shemi & Piran, 1990; 
Paczyhski, 1990). The kinetic energy can be recovered as radiation only if these 
particles are decelerated by the ISM (Meszaros, & Rees, 1992) or by internal 
shocks (Rees & Meszaros 1994, Narayan, Paczyhski & Piran, 1992). In retro- 
spect, once this is realized, one can imagine GRB models in which the fireball 
is replaced by another (unknown) non-thermal acceleration mechanism but the 
radiation is still emitted due to slowing down of the ultra-relativistic particles. 
It is worthwhile, therefore, to explore the nature of the interaction between the 
ultra-relativistic particles and the ISM. We show that a careful analysis of this 
interaction changes some of the previous results and it may shed a new light on 
the expected temporal structure in GRBs. 

We examine, first, in section 2 the planar shock problem. Spherical effects 
play, however, a crucial role in the realistic situation and we consider them in 
section 3. We discuss the observational implications to GRBs in section 4. 

2. Planar Symmetry 

Consider a slab of ultra-relativistic cold dense matter with a Lorentz factor 7 ^> 
1 that hits a stationary cold interstellar medium (ISM). Two shocks form: a 
reverse shock that propagates into the dense relativistic shell, reducing its speed 
and increasing its internal energy, and a forward shock that propagates into the 
Interstellar medium giving it relativistic velocities and internal energy. A contact 
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discontinuity separates the shocked shell material and the shocked ISM material 
(Meszaros & Rees 1992; Katz, 1994). 

There are four regions in this system: the ISM (1), the shocked ISM (2), 
the shocked shell material (3) and the unshocked shell material (4). The ISM is 
at rest relative to the observer. Velocities and their corresponding Lorentz 
factors 7; = (1 — fif) , distances and time are measured relative to this frame. 
Thermodynamic quantities: pi and e 4 - (particle number density, pressure and 
internal energy density) are measured in the fluids' rest frames. The ISM and 
the unshocked shell are cold and therefore: e\ = e 4 = 0. The shocked material is 
extremely hot and therefore: p 2 = e 2 /3 and p 3 = e 3 /3. 

For 7 = 74 ^> 1 the equations governing the shocks are (Blandford & McKee, 
1976): 

e 2 /n 2 m p c 2 = 72 - 1 = 72 ; n 2 /n 1 = A^ 2 + 3 = 472 , (1) 

e 3 /n 3 m p c 2 = 73 - 1 ; n 3 /n 4 = 473 + 3 , (2) 

where m p is the protons rest mass. The approximations in equation 1 used only 
the fact that 74 ^> 1 and therefore 72 3> 1. No assumption was made about 73, 
the Lorentz factor of the motion of the shocked material in region 3 relative to 
the unshocked shell in region 4. 

Equality of pressures and velocities along the contact discontinuity yields: 

e 2 = e 3 ; 73 = (74/72 + 72/7*0/2 (3) 

The solution for 72 depends only on two parameters 7 and / = n±jri\. The energy, 
pressure and density also depend linearly on a third parameter, the external 
density ri\. A forth parameter, A, the width (in the observer's frame) of the 
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ultra-relativistic shell determines the time it takes the reverse shock to cross the 
shell, t&: 

A / 7n 4 \ 

There are two simple limits of equations 1-4 in which the reverse shock is either 
Newtonian or ultra-relativistic (the forward shock is always ultra-relativistic if 
7 3> 1 and / > I/7 2 ). If 7 2 3> / the reverse shock is ultra-relativistic (73 ^> 1) 
and: 

y/2 7 1/2 / 1/4 

73 = 72f^ ; 72 = 73 = ^~ (5) 

In this case almost all of the initial kinetic energy is converted by the shocks into 
internal energy (73 <C 7)- Therefore the process is over after a single passage of 
the reverse shock through the shell. The relevant time scale for energy extraction 
is the shell crossing time: 

t A = A 7 /f/2c . (6) 

The internal energy densities in the shocked shell and in the shocked ISM is the 
same (see equation 3) and since both shocked regions have comparable width 
they release comparable amounts of energy. The ISM mass swept by the forward 
shock at the time that the reverse shock crosses the shell is ~ J -1 / 2 of the shell's 
mass. This is larger than the simple estimate given by Meszaros & Rees (1992) 
of ~ 7 _1 . At an earlier time when a mass of 7 _1 was swept, the reverse shock 
interacted only with a small fraction of the shell (~ v7/7 ^ 1) an d most of the 
energy was still the kinetic energy of the unshocked shell. 

If / 3> 7 2 the reverse shock is Newtonian (73 - 1 < 1) and: 

73-1^ 4 7 2 r 1 /7 = 2 e <l ; 72 = 73 = 7(1-^). (7) 
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The shock crosses the shell at: 



yfi/UAjy/f/c, 




which is surprisingly similar (up to a constant) to the ultra-relativistic limit 
expression. 

The reverse shock converts only a fraction 7/v7 <C 1 of the kinetic energy 
into internal energy. It is too weak to slow down the shell effectively and most 
of the initial energy is still kinetic energy when this shock reaches the inner edge 
of the shell. At this stage a rarefaction wave begins to propagate towards the 



contact discontinuity. This wave propagates at the speed of sound A/4p 3 /3n 3 m. 



the same order of magnitude as the shock crossing time t&. It is then reflected 
from the contact discontinuity and a second, weaker, shock wave forms. A quasi- 
steady state slowing down solution forms after a few crossings like this (Sari & 
Piran, 1995). Using momentum conservation, the total slowing down time can 
be estimated by ~ A^n^m-pc/ p 2 ~ Af /cy. During this time the forward shock 
collects a fraction ~ 7 _1 of the shell's rest mass, which is the same as the original 
estimate of Meszaros & Rees (1992). In contrary to the relativistic case, there 
are two time scales now: the shock (or rarefaction) crossing time, and the 
total slowing down time. 

In the realistic situation the ISM density is probably inhomogeneous. Consider 
a density jump by a factor /' over a distance Iism- The forward shock propagates 
into the ISM with a density n\ as before and when it reaches the position where 
the ISM density is n\f a new shock wave is reflected. The solution of this 
problem requires the application of equations similar to equations 1-4. This 




and it reaches the contact 




which is of 
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shock is reflected again of the shell. Similar analyses shows that the reflections 
time is ~ lisml^c^/J 1 and after these reflections the pressure and time scales are 
as if the ISM was homogeneous with a density n\f (Sari & Piran, 1995). 

Finally, we mention the possibility of internal shocks inside the shell (Rees 
& Meszaros 1994; Narayan, Paczyhski & Piran, 1992). These may form when 
faster material overtakes slower material. If the Lorentz factor varies by a factor 
of ~ 2 over a length scale 6R < A then the time for these shock to from is 
~ SR~/ 2 /c < A7 2 / c. This time scale is shorter than the slowing-down time 
scale and therefore internal shocks appear before considerable deceleration in the 
Newtonian case. In the relativistic case considerable deceleration occurs before 
internal shocks unless 6R <C A. 

3. Spherical Considerations 

The main difference between spherical and planar symmetries is that in a spherical 
system the density ratio / = n 4 /rii decreases with time. Initially f /~f 2 3> 1 and 
the reverse shock is Newtonian. The energy conversion depends critically on 
the question whether this shock become relativistic before the kinetic energy is 
extracted from the shell. This depends on the ratio of two radii: where 
/ 1^ 2 = 1 and the reverse shock becomes relativistic and where the reserve 
shock crosses the shell. Two other important radii are: i? 7 where the forward 
shock sweeps a mass Af/7 (M is the shell's rest mass) and R s where the shell 
begins to spread if the initial Lorentz factor varies by order 7 (Meszaros, Laguna, 
& Rees, 1993; Piran, Shemi & Narayan, 1993; Piran, 1994). Note that R s is also 
an upper limit for the the location of internal shocks since 6R < A. 
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There are two intrinsic length scales in this problem. The first is A, the width 
of the relativistic shell. The second is the Sedov length, / = (E /nirripC 2 ) 1 / 3 , which 
is familiar from SNR theory. The ISM rest mass within / 3 equals E/c 2 . Typical 
GRB parameters yield, / 10 18 cm, which is similar to SNR value. We estimate 
the mass of the shell using, M = E j^c 2 , and obtain i? 7 = l/^f 2 / 3 . For a shell that 
propagates with a constant width / oc R~ 2 and we can express the other radii 
(omitting here and in the rest of the discussion factors of order of unity) in terms 
of /, A and 7: R N = I 3 / 2 / A 1 / 2 ^ 2 ; R A = \ 3 > A A 1 / 4 and R s = A7 2 . Conveniently, 
the four critical radii are related by one dimensionless quantity: 

C = (//A) 1/2 7 - 4 / 3 . (9) 

and 

R N /C = R 1 = ^R A = C 2 R S ■ (10) 

Two possibilities exist: 

1. £ > 1 - the Newtonian case: R s < R A < i? 7 < and shock reaches the 
inner edge of the shell while it is still Newtonian. Most of the energy is extracted 
during a steady state deceleration phase described in the previous section. Since 
R s is smaller than all other radii spreading migth be important. If the shell is 
spreading then A in the above expressions should be replaced by -R/7 2 . This 
delays the time at which the reverse shock reaches the shell and decreases the 
shell's density. These effects lead to a triple coincidence: = i? 7 = with 
£ ~ 1 and a mildly relativistic reverse shock during the period of effective en- 
ergy extraction. Without spreading only a small fraction of the total energy is 
converted to thermal energy in the reverse shock. With spreading both shocks 
convert comparable amounts of energy. 
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2. £ < 1 - the relativistic case: < R'J < Ra < R s - The reverse shock 
becomes relativistic before it crosses the shell. Only a small fraction of the energy 
is converted at i? 7 and the kinetic energy is converted into internal energy only 
at i?A- R s is larger than all other radii and spreading is unimportant. It is 
interesting to note that in this limit ^2{Ra) ~ (//A) 3 / 8 is independent of the 
initial Lorentz factor 7 and it is only weakly dependent on other parameters. 
This might have an important role in the fact that the observed radiation always 
appears as low energy 7-rays. 

Neither the internal shocks nor the ISM inhomogeneity time scales are affected 
by these spherical considerations. The former depends only upon 6R } A and 7 
and the latter depends only on Iism and 7. Both are constant throughout the 
spherical expansions. 

4. Observational Implications to GRB 

We estimate now the relevant parameters for GRBs and examine the possible 
relation between the observed time scales and the hydrodynamic time scales. We 
assume that the shocked material (either region 2 or 3) emits the radiation on 
a time scale shorter than the hydrodynamic time scales. A simple estimate of 
synchrotron cooling rate (assuming equipartition of the magnetic field energy) is 
consistent with this assumption. 

The total energy of the bursts can be estimated directly from the observed 
fluxes (assuming cosmological distances) as: E = I0 51 ergs. The ISM density has 
a typical value of: n\ = Iparticle/cm 3 . Only the ratio of these two quantities 
appears in our considerations and it determines the Sedov length / I0 18 cm. 
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The values of A and 7 are more ambiguous. It is known that 7 > 100 in order for 
the shell to be transparent for 7-rays (Fenimore, Epstein & Ho, 1993; Woods & 
Loeb, 1995; Piran, 1995). A similar constraint can be obtained from the observed 
duration of the bursts and we adopt 7 = 10 3 as our canonical value. The width 
of the shell is highly uncertain and relativistic effects allow it to be several order 
of magnitude larger than the common canonical value A = 10 7 cm. For these 
canonical parameters £ = 30 > 1, corresponding to a Newtonian reverse shock. 
Nevertheless a value of £ = 0.1 < 1 is also possible with reasonable parameters 
(for example A = 10 9 cm and 7 = 10 4 ). Therefore both relativistic and Newtonian 
reverse shock are possible. 

The bursts' duration is determined by the slowing down time of the shell. 
The emitting region moves towards the observer with a Lorentz factor 72. Two 
photons that are emitted with a time delay dt will be detected with a time delay 
dt/^- Additionally, an observer detects radiation from a region with an angular 
size 72 -1 - A photon emerging from an angle away from the center of a region 
with a radius R will be detected at a time Rj^\c after a photon that emerges 
from the center (Katz, 1994). Thus, given a typical radius of energy conversion, 
R e the observed time scale is: 

{A/c if £ < 1 (Relativistic); 

(H) 
i? 7 / 7 2 c ~ // 7 8 / 3 c if £ > 1 (Newtonian) . 

This time scale ranges from ~ 1msec, for A = 3 X 10 7 cm and 7 = 10 4 , to 

~ lOOsec for 7 = 10 2 and A = 10 13 cm (see fig. 1). In the Newtonian regime 

(£ > 1) the observed time scale depends only on 7 while it depends only on 

A in the relativistic (£ < 1) case. The observed durations of the brightest 30 

bursts limit 7 to 100 < 7 < 10 4 with a typical value of 500 and it limits A to 
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A < 3 x 10 12 cm. 

In the Newtonian case there appears a second time scale: the crossing time 
ol the shell by the reverse shock. The corresponding observed time scale: 

t~/ 3 / 4 A 1 / 4 / 7 2 c = t o6s /^, (12) 

is comparable to the measured time scale of variability in the bursts (10msec- 
lOsec). However, it not clear if this time scale has any observational implications 
since only a small fraction of the energy is emitted by the reverse shock in this 
case (it might though produce primary photons that will be Compton scattered 
later by the forward shock). Furthermore, spreading prolongs this scale so that 

t ~ t bs- 

Another, more likely, source of the variability is inhomogeneity in the ISM. If 
the length scale of the inhomogeneity is Iism and the density varies by one order 
of magnitude then the time scale for the observed variability will be: 

tvar ~ hsml^llc (13) 

This time scale can be sufficiently short to produce the observed variability if 
hsM is sufficiently small. 

About 3% of the bursts contain precursors: weaker bursts that proceed the 
main burst. Two additional time scales appears here: the precursor's duration 
At pre , and its separation from the main burst At pre _ ma8m . A natural explanation 
for the precursor phenomenon, within this model, is that it arises from internal 
shocks that take place at R 6R^ 2 < A7 2 while the main burst originates from 
the interaction with the ISM. (Meszaros, & Rees; 1994) proposed that internal 
shock produce the main GRB while the interaction with the ISM produces the 
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delayed GeV photons observed in some bursts). The duration of the precursor is 
At pre = A/c, which requires values of A as high as 10 10 — 10 12 cm to produce the 
observed precursors of 1 — 100 sec. If £ > 1 the main bursts will have a duration 
At f, s //7 8/,3 c = <^ 2 A. This will also be the typical separation At pre _ ma8m 
between the precursor and the main burst. We expect a time delay between the 
precursor and the main burst which will be comparable to the duration of the 
main burst. Note that a correlation of the form At pre _ ma8m 4.5At 6 s exists (but 
was not reported) in the data of Koshut et. al., (1995). If £ < 1 then precursors 
do not occur (unless 6R <C A) This is in agreement with the lack of observed 
precursors in short bursts. 

5. Conclusions 

We have calculated the hydrodynamic time scales of shocks during the interac- 
tion between an ultra-relativistic shell and the ISM. These time scales depend on 
the shock conditions which in turn depend only on energy and momentum con- 
servations. Hence, we believe that these time scales are robust and independent 
of the unknown microphysics that takes place in these shocks. We find that with 
reasonable astrophysical parameters these time scales are in a good agreement 
with the observed time scales in GRBs. Our analysis shows that there are two 
kinds of shocks: Newtonian and Relativistic. The difference between them might 
correspond to the difference between the observed short and long bursts (Kou- 
veliotou et. al. 1993). Finally, we suggest that precursors might be emitted due 
to internal shocks within the ultra-relativistic shell while the main burst emerges 
later from the interaction with the ISM. 
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Figure Captions 

Figure 1: The observed duration of the burst (dashed curve in the Newtonain 
regime and dotted curve in the relativistic regime) as function of 7 and A for 
/ = 10 18 cm (E = 10 51 erg and n\ = lcm~ 3 ). The thick solid line (£ = 1) separates 
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the Newtonian (lower left) and the Relativistic (upper right) regions. Spreading 
drives all Newtonian cases to the £ ~ 1 line. 
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